Universal thermodynamic behavior of supercritical matter 
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Since their discovery in 1822, the supercritical fluids have been of enduring interest, and have 
started to be deployed in a number of important industrial applications. Theoretical understanding 
of the supercritical state is lacking, and is seen to limit further industrial deployment. Here, we 
study thermodynamic properties of the supercritical state, and discover that specific heat shows 
a crossover between two different regimes, an unexpected result in view of currently perceived 
homogeneity of supercritical state in terms of physical properties. We subsequently formulate a 
theory of system thermodynamics above the crossover and find good agreement between calculated 
and experimental specific heat with no free fitting parameters. We derive a power law and analyze 
supercritical scaling exponents in the system. In this theory, energy and heat capacity are governed 
by the minimal length of the longitudinal mode in the system only, and do not explicitly depend on 
system-specific structure and interactions. In this sense, supercritical properties exhibit universality 
not hitherto anticipated. 
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PACS numbers: 05.70.Fh, 65.20.Jk, 89.75.Da 

Statistical mechanics is the art of predicting the be- 
havior of a system with a large number of degrees of 
freedom, given the laws governing its microscopic behav- 
ior. The statistical description of liquids, in comparison 
with the solid and gas phases, is incomplete. The prob- 
lem of formulating a rigorous mathematical description of 
liquids has always been regarded as much more difficult 
than that of the kinetic theory of gases or solid-state the- 
ory, stimulating the ongoing research (ll-ITo|. Due to the 
simultaneous presence of strong interactions and large 
atomic displacements, common models and approxima- 
tions used for gases and solids, do not apply to liquids. 
For this reason, liquids do not generally fall into any sim- 
ple classification, and have been mostly treated as general 
many-body systems as a result. 

In recent years, a significant effort has been devoted 
to investigation of various properties of supercritical flu- 
ids [ll| - [l5j . This has been an exciting field with a long 
history since 1822 when Baron Charles Cagniard de la 
Tour discovered supercritical fluids while conducting ex- 
periments involving the discontinuities of the sound in 
a sealed cannon barrel filled with various fluids at high 
temperature [3] • More recently, supercritical fluids have 
started to be deployed in several important applications, 
ranging from the extraction of floral fragrance from flow- 
ers to applications in food science such as creating decaf- 
feinated coffee, functional food ingredients, pharmaceu- 
ticals, cosmetics, polymers, powders, bio- and functional 
materials, nano-systems, natural products, biotechnol- 
ogy, fossil and bio-fuels, microelectronics, energy and en- 
vironment 14, ]j| 171. Much of the excitement and in- 



terest of the past decade is due to the enormous progress 
made in increasing the power of relevant experimental 
tools [18l - l2l| . The development of new experimental 



methods and improvement of existing ones continues to 



play a important role in this field [22 - 27 1 , with recent re- 
search focusing on dynamic properties of fluids 

High density and high thermal motion are two main 
properties responsible for efficient cleaning, dissolving 
and extracting abilities of supercritical fluids in the above 
industrial applications. From the point of view of prac- 
tical applications, supercritical fluids have got the best 
of both worlds: high density comparable to ordinary liq- 
uids and solids and high thermal motion and diffusiv- 
ity approaching that of gases. Notably, it is this very 
combination that presents a formidable problem to the 
theory: high density and strong interactions mean that 
theories and approximations used for dilute gases do not 
apply [34j]. Early Enskog and related kinetic approaches 
to gases were followed by more extensive developments, 
yet they do not adequately describe dense systems with 
strong interactions and many-body correlations [3| such 
as supercritical fluids. One general issue with extending 
gas-like approaches to fluids was noted earlier: in a sys- 
tem with strong interactions, the system energy strongly 
depends on the type of interactions, and is therefore 
system-specific, ruling out the possibility to develop a 
theory universally applicable to many fluids, in contrast 
to gases and solids [351 ] . 

In addition to theoretical challenges, the lack of fun- 
damental understanding is seen as an obstacle towards 
wider deployment of supercritical fluids in industrial ap- 
plications, primarily due to the absence of guidance re- 
garding pressure and temperature at which the desired 
properties are optimized as well as the possibility to use 
new systems |14j . 

In this paper, we focus on the thermodynamic proper- 
ties of the supercritical state. On the basis of molecular 
dynamics simulations, we find that specific heat shows a 
crossover between two different dynamic regimes of the 
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low-temperature rigid liquid and high-temperature non- 
rigid gas-like fluid. The crossover challenges the currently 
held belief that no difference can be made between a gas 
and a liquid above the critical point and that the super- 
critical state is homogeneous in terms of physical proper- 
ties [III. We subsequently formulate a theory of system 
thermodynamics and heat capacity above the crossover. 
In this theory, energy and heat capacity are governed 
by the minimal length of the longitudinal mode in the 
system only, and do not depend on system-specific struc- 
ture and interactions. We further study the predicted 
relationship between supercritical exponents of heat ca- 
pacity and viscosity. A good agreement is demonstrated 
between calculated and experimental data for noble and 
molecular supercritical fluids with no free fitting param- 
eters. 



RESULTS 
Dynamic crossover of the specific heat 

We start with molecular dynamics simulations of a 
model liquid. Our primary aim here is to show that spe- 
cific heat, cy, shows a crossover in the supercritical region 
of the phase diagram. This result is unexpected in view 
of currently perceived homogeneity of supercritical state 
in terms of physical properties. 

Using molecular dynamics simulations [37l |. we have 



simulated the binary Lennard- Jones (LJ) fluid [38|. We 
have simulated the system with 64000 atoms using 
constant- volume (nve) ensemble in the wide temperature 
range (see FigfT]) well extending into the supercritical re- 
gion (the critical temperature of this system is T c ~ 160 
K). From the energy of the system E at each tempera- 
ture, we calculate constant-volume specific heat, cy, as 

CV = jf§ (ka = 1). 

We observe that cy decreases steeply from the solid- 
state value of about 3fce at low temperature to approxi- 
mately 2&B around 2000 K. The steep decrease is followed 
by crossing over to a considerably weaker temperature de- 
pendence. This crossover is a new effect not reported in 
previous MD simulations. We further observe that the 
crossover takes place around cy « 2fce. This value of 
cy, cy = 2fce is non-coincidental, and corresponds to the 
crossover taking place across the Frenkel line 39l 4pj] as 
discussed below. 

Crossing the Frenkel line corresponds to the qualita- 
tive change of atomic dynamics in a liquid. In liquids, 
atomic motion has two components: a solid-like, quasi- 
harmonic vibrational motion about equilibrium locations 
and diffusive gas-like jumps between neighboring equi- 
librium positions. As the temperature increases, a par- 
ticle spends less time vibrating and more time diffusing. 
Eventually, the solid-like oscillating component of motion 
disappears; all that remains is the gas-like ballistic mo- 
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FIG. 1. Calculated cv showing the crossover and continuous 
dynamical transition around cv ~ 2, (fee = 1). The crossover 
takes place between different dynamical regimes of the rigid 
liquid and non-rigid supercritical fluid. System density is 2072 
kg/m 3 , approximately four times larger than critical density 
of liquid Ar. 



tion. That disappearance, a qualitative change in particle 
dynamics, corresponds to crossing the Frenkel line, the 
transition of the substance from the liquid dynamics to 
the gas dynamics. This transition takes place when liquid 
relaxation time t (t is liquid relaxation time, the average 
time between consecutive atomic jumps at one point in 
space [4l| ) approaches its minimal value, rp, the Debye 
vibration period. As recently discussed 39|, 4(| , crossing 
the Frenkel line is accompanied by qualitative changes 
of most important properties of the system, including 
diffusion, viscosity, thermal conductivity and dispersion 
curves [39L l40l| . One of these properties is the crossover of 
the speed of sound and appearance of the high-frequency 
"fast" sound 2l|, |42j-|44( . At the microscopic level, the 



Frenkel line can be identified by the loss of oscillations in 
the velocity autocorrelation function fiol ] , a point studied 
in detail elsewhere. 

Fig. 1 implies that in addition to the dynamic proper- 
ties discussed above, the thermodynamics of the system 
changes at the Frenkel line too. This is an important 
general insight regarding the behavior of the supercriti- 
cal state. 

The initial steep decrease of cy from about 3fce to 2fce 
can be quantitatively explained by the progressive loss 
of two transverse waves with frequency w > 1/r (SSl. 
Physically, this picture is based on Frenkel's prediction 
that on time scale shorter than t, liquid is a solid, and 
therefore supports rigidity and solid-like transverse waves 
at short times, or at frequency larger than u) > 1/r [48| . 
When t approaches its maximal value, td, the liquid 
can not sustain transverse waves at any frequency [49j |. 
Consequently, the potential energy of the system is due 
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to the longitudinal mode only, giving the total energy 



of 2 AT and specific heat of 2fc B [45|, |47[. Hence, the 
decrease of cy from about 3/cb to 2k B corresponds to 
the region of a "rigid" liquid where short-time solid- like 
rigidity and high-frequency transverse waves exist. On 
the other hand, the liquid is unable to sustain transverse 
waves at any available frequency above the Frenkel line. 
Instead, the liquid enters a new dynamic "non-rigid" gas- 
like regime where oscillatory component of particles is 
lost and the motion becomes purely collisional as in a 
gas. 

We therefore need to develop a new theory capable of 
describing thermodynamics of supercritical matter above 
the Frenkel line where the system enters the new dynamic 
regime and where cy falls below 2k B and approaches the 
ideal-gas value of 3/2fcB at high temperature (see Fig. 
®. 



Thermodynamic theory of supercritical state above 
the Frenkel line 



We now focus on the theory of the non-rigid gas-like 
liquid above the Frenkel line, and add a new proposal 
regarding how the system energy can be evaluated. As 
temperature rises in the ballistic gas-like regime and ki- 
netic energy increases, the mean free path I, the average 
distance between particle collisions, increases. At the 
Frenkel line where the ballistic regime starts, I is com- 
parable to interatomic separation. In the limit of high 
temperature where the particle's kinetic energy is much 
larger than potential energy, I tends to infinity as in the 
non-interacting ideal gas. Our proposal is that I deter- 
mines the shortest wavelength of the longitudinal mode 
that exists in the system, A, because below this length 
the motion is purely ballistic and therefore can not be 
oscillatory: X — I. On the other hand, the longitudinal 
modes with larger wavelength are supported, and repre- 
sent the excitations existing in the supercritical system. 

We note that the existence of long- wavelength longitu- 
dinal waves in a gas, sound, is well-known. What is new 
here is that we propose that the contribution of the lon- 
gitudinal waves to the energy of the gas-like supercritical 
system starts from very short wavelengths comparable 
with interatomic separation a. In this sense, we are ex- 
tending the solid state concepts (e.g. short-wavelength 
solid-like phonons with Debye density of states, see be- 
low) to the new area of gas-like supercritical state where 
these ideas have not been hitherto contemplated. We also 
note that experimentally, it is well established that dy- 
namics in sub-critical liquids shows solid-like character, 
in that liquids can sustain high-frequency propagating 
modes down to wavelengths on the atomic scale, with 
solid- like dispersion relations (49l |. 

Here and elsewhere, our discussion of liquid vibrational 
states includes an important point. Namely, these vi- 



brational states are not harmonic plane waves that de- 
cay in disordered media. Instead, we refer to liquid vi- 
brational modes as non-decaying solutions, eigenstates, 
of the secular equation involving the force matrix con- 
structed from the amorphous liquid structure. Experi- 
mentally, the dispersion relation for these waves is quasi- 
linear as for plane waves, leading to the quadratic density 
of states g(u>) oc w 2 . A detailed discussion of this point 
is forthcoming. 

In the proposed theory, the energy of the non-rigid 
supercritical fluid per particle includes the contribution 
from the kinetic energy, K — |fc B T and the potential en- 
ergy of the longitudinal phonons with wavelengths larger 
than I. Using the equipartition theorem (Pi) = ^rf^, 
where (E{) is the energy of the longitudinal phonons, we 
write 



E = ^k B T + i J e(u, T)g(u)du + E anh 



(1) 



where the upper integration limit ljq is given by the short- 
est wavelength in the system, A: ojq = ^-c, c is the speed 
of sound, e(u!,T) is the mean energy of harmonic oscil- 
lator, e(u,T) = ?f + e J^_ i ; or e(T) = k B T in the 
classical case, and E an h is the anharmonic contribution 
to the phonon energy. 

The second term in Eq. ((SJ can be calculated using the 
Debye density of states, g(w) = -4-w 2 , giving k B T^-, 

3 

or ksTjs, where a is the interatomic separation. The 
use of the quadratic density of states is supported by the 
experimental evidence showing solid-like quasi-linear dis- 
persion relationships in supercritical fluids 21 , H3| similar 
to the subcritical liquids. E an h can be evaluated in the 
Griineisen approximation from the softening of phonon 
frequencies with temperature, with the result that the 
energy is modified as E — > E (l + ^) where a is the 
coefficient of thermal expansion [46L T47j| . Then, the en- 
ergy of non-rigid gas-like supercritical fluid becomes: 



E=-k B T- 



1 



-aT 



1 a 3 
2 te7 > 



(2) 



We observe that when A s» a at the Frenkel line, Eq. 
([2]) gives E = 2k B T and the specific heat of 2 (here we 
neglect the small term cx aT). This corresponds to the 
crossover of cy in Figure 1 to the gas-like regime as dis- 
cussed above. When A increases on temperature increase 
in the gas-like regime, Eq. ^ predicts that cy tends to 



§, the ideal-gas value as expected. 



Comparison with experimental data 

As discussed above, A is given by the particle mean 
free path in the non-rigid gas-like regime, and can there- 
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fore be calculated from n = ^puX, where n is viscosity, p 
is density and u is average velocity. Therefore, Eq. ([2]) 
provides an important relationship between the energy of 
the non-rigid supercritical liquid and its gas-like viscos- 
ity. We now check this relationship experimentally, by 
comparing the specific heat cy — predicted by Eq. 
([2]) and the experimental cy. 

We have used the National Institute of Standards and 
Technology (NIST) database [HI . We aimed to check our 
theoretical predictions and selected the isochoric data of 
several supercritical noble and molecular liquids in a wide 
range of temperature and density. We note that A in Eq. 
([2]) changes in response to both temperature and den- 
sity: A increases with temperature and decreases with 
density. Practically, the range of isochoric data in the 
NIST database is fairly narrow in terms of density com- 
pared to the range of temperature. We therefore analyze 
the temperature behavior of c v and rj along several iso- 
chores. Our choice of liquids is dictated by the availabil- 
ity of isochoric data in the supercritical region. For each 
density, we fit experimental viscosity to r\ = Aq + A\T A2 , 
calculate A from n = | pu\ and subsequently use A in Eq. 
([2]) to calculate cy. 

We note that Debye model is not a good approxima- 
tion in molecular liquids where the frequency of intra- 
molecular vibrations considerably exceeds the rest of fre- 
quencies in the system (3340 K in N2 and 3572 K in CO). 
However, the intra-molecular modes are not excited in 
the temperature range of experimental cy (see Figures 
34). Therefore, the contribution of intra-molecular mo- 
tion to cv is purely rotational, c rot . On the other hand, 
the rotational motion is excited in the considered tem- 
perature range, and is therefore classical, giving c rot — R 
for linear molecules in N2 and CO. Consequently, cy for 
molecular liquids shown in Figure 4 correspond to heat 
capacities per molecule, with c rot subtracted from the 
experimental data. 

We also note that experimental isochoric cy is affected 
by lambda-like critical anomalies (see Figs. 2, 3) be- 
cause the isochoric NIST data do not extend far above 
the critical point. Here, we do not consider critical effects 
related to phase transitions, and therefore fit the data at 
temperatures that are high enough to be affected by the 
lambda-anomaly at the phase transition. In Figures 2 
and 3 we observe good agreement between experimental 
and predicted cy, in view of (a) 2-5% uncertainty of ex- 
perimental cy and 77 [HH , (b) approximations introduced 
by the Debye model and (c) increased curvature of cy at 
low temperature due to proximity of lambda-anomalies 
that are not taken into account by the theory. Notably, 
the agreement is achieved without using free fitting pa- 
rameters because p, a, u and a are fixed by system prop- 
erties. Values of these parameters used in Eq.([2]) are in 
good agreement with their experimental values. 

Another way to compare our theory and experimental 
data is to study supercritical exponents of specific heat 
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FIG. 2. Experimental and calculated cv {ks = 1) for non- 
rigid supercritical fluid Neon. Experimental c v and 77 are 
taken from the NIST database (51| at density 700 kg/m 3 . 
Value of a used in the calculation is 3.5 x 10~ 3 K _1 . 



and viscosity. Indeed, if n in the non-rigid gas-like su- 
percritical region can be approximated as a power law, 
ij oc T 7 , then Eq. © makes two predictions. First, E 
and cy should also follow power laws. Second, Eq. ([2]) 
provides a specific relationship between the scaling expo- 
nents of n on one hand and E and cy on the other hand, 
the relationship that we check below. 

Apart from comparing theory and experiments, study- 
ing the supercritical scaling exponents and discussing 
their possible universality is interesting in the wider con- 
text of scaling behavior of physical properties. In the 
area of phase transitions, the scaling behavior idea has 
been a crucial element in the subject of liquids and other 
systems from the nineteenth century onwards JEsLCriti- 
cal points occur in a great variety of systems [53j-|56j , yet 
there is a considerable degree of similarity (universality) 
in the way in which systems approach the critical point 
[57j ]. Here, the calculation of critical exponents serves as 
one of the main test of the theories. 

Once analyzed within our theory below, experimen- 
tal temperature dependence of n, u, E and cy [5l| can 
be fairly well approximated by the power law above the 
Frenkel line. Therefore, we use the following power-law 
relationships: 



rj = const x T 7 
1 

u = const x T 2 

E= 1.5 X T + const x T s 
const 

Cy = 1.5 



(3) 
(4) 
(5) 

(6) 



From rj = \pu\ (k B = 1), A = const x T 7 ° 5 . Using 
A in Eq. (J5|) gives for the energy (neglecting the small 
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FIG. 3. Experimental and calculated cv (ks = 1) for noble non-rigid supercritical fluids. Experimental c v and r\ are taken 
from the NIST database [5l|] at different densities as shown in the Figure. Values of a used in the calculation are 3.5 x 10~ 3 
K~ l (Ne), 3.3 x 10" 3 K' 1 (At), 1.8 x 10" 3 K" 1 (Kr) and 8.2 x 1CT 4 KT 1 (Xe). The uncertainty of both experimental heat 
capacities and viscosities is about 2-5% [5l|]. Insets show viscosity fits. 
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FIG. 4. Experimental and calculated cv (&b = 1) for molecular non-rigid supercritical fluids. Experimental c v and r\ are taken 
from the NIST database [5l[. Values of a used in the calculation are 6.5 x 10" 3 K _1 (CO) and 8.5 x 10" 3 K _1 (N 2 ). The 
uncertainty of both experimental heat capacities and viscosities is about 2-5% [HJ. Insets show viscosity fits. 



oc aT term): 
and specific heat: 



Supercritical Fluid < £ > < jth > < ~fex P > 



E oc T 2 - 5 " 37 



c v oc (2.5 - 37)T 



1.5-37 



(7) 



(8) 



Then, from Eqs. ©, © and © we find the following 
relationships between the power-law scaling exponents: 



2.5 -3 7 = 8 



37 - 1-5 = i 



(9) 



(10) 



We note that Eqs. © and flSJ, together with the gen- 
eral requirement for cy to be positive and the experi- 
mental requirement for cy to decrease with temperature, 
imply that 7 should be in the range \ < 7 < |. This is 
the case for the experimental data we examined, as will 
be shown below. We also note that 7 = \ corresponds 
to the non-interacting ideal gas and r\ oc T 0,5 . 

We now check the Eq. (ITU|) experimentally. Using the 
NIST database [HT| and Eq.©, we calculate £ from the 

slope of log (cy — 1.5) vs. log ( T ^ Tp ) , where Tp corre- 
sponds to the temperature on the Frenkel line and 1.5 is 
the asymptotic gas value of specific heat. Using Eq. (ITUl) . 
we calculate the predicted theoretical value jth and com- 
pare it with the experimental 7 ea; p obtained from fitting 
the experimental viscosity to 77 oc T 7ex p. 

In Table U we show £, ~f exp and 7^, averaged over sev- 
eral data sets taken along the isochores at several dif- 
ferent densities. We observe the overall good agreement 
between j exp and jth- We further observe that the power 
exponent of specific heat, £, is close to 0.2 for different 
systems. In our theory, the near universality of temper- 
ature behavior of cy is due to temperature dependence 
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TABLE I. Power-law exponents £, y e x P and 7th, averaged over 
different densities: Ne(600 kg/m 3 , 700 kg/m 3 , 900 kg/m 3 ), 
Ar (700 kg/m 3 , 800 kg/m 3 , 900 kg/m 3 ), Kr(1100 kg/m 3 , 
1300 kg/m 3 , 1500 kg/m 3 ), Xe(1300 kg/m 3 , 1500 kg/m 3 , 1700 
kg/m 3 ), CO (350 kg/m 3 ) and N 2 (320 kg/m 3 ). 



of A. This point is discussed in the next section in more 
detail. 



DISCUSSION 

Our first important observation in this work is that 
contrary to the current belief, the thermodynamic prop- 
erties of the supercritical state are not homogeneous. In- 
stead, the specific heat shows a crossover related to the 
change of particle dynamics, which we attributed to the 
recently introduced Frenkel line. 

We have subsequently focused on thermodynamic 
properties of supercritical fluids above the Frenkel line. 
Here, we faced the problem of strong interactions, the 
long-persisting challenge in condensed matter physics. 
Indeed, strong interactions imply that approximations 
used for dilute gases do not apply to real dense liquids 
34|. If we consider realistic strong interactions (assum- 
ing that interactions are known and can be represented 
in analytical form, the assumption that is valid for a rel- 
atively small number of simple systems only) and struc- 



7 



tural correlations that often include those beyond two- 
body correlations, we quickly find that the problem be- 
comes intractable. Further, strong interactions, coupled 
with their specificity in different systems, have been sug- 
gested to preclude the calculation of en ergy and heat ca- 
pacity in general form from the outset [35[. 

In this paper, we addressed the problem in a different 
way, by substituting all potentially complicated effects 
of interactions and structural correlations by one physi- 
cal quantity, the minimal wavelength of the longitudinal 
mode in the system A. This has enabled us to rational- 
ize the experimental behavior of cy as well as provide the 
relationship between different physical properties and ex- 
perimental outcomes (e.g. relationship between cy and 
77). Notably, our approach unveils universality of ther- 
modynamics of supercritical state in the following sense. 
First, cy does not explicitly depend on system details 
such as structure and interactions but on A only. Fluids 
may have very different structure and interactions, yet 
our theory predicts the similarity of their thermodynamic 
behavior as long as A behaves similarly in those systems. 
Second, and more specifically, our approach predicts that 
supercritical scaling of thermodynamic properties such 
as heat capacity is governed by viscosity scaling. Con- 
sequently, similar temperature scaling of viscosity gives 
universal temperature scaling of thermodynamic proper- 
ties. 
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